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, After the submission of the paper, three strong earthquakes with magnitude around 6.0-units 

occurred on October 17 and October 20, 2005, with epicenters in the Aegean Sea, at a distance only 
100km from MYT station at which the intense signals Mi to M4 -analyzed in the main text- have 
been recorded. This confirms experimentally the proposed criterion we used for the classification of 
these signals as Seismic Electric Signals (SES). Moreover, we show that, if we follow the procedure 
described in [P.A. Varotsos, N. V. Sarlis, H. K. Tanaka and E. S. Skordas Phys. Rev. E 72, 
041103 (2005)], the analysis in the natural time of the seismicity after the SES initiation allows the 
estimation of the time window of the impending earthquakes with very good accuracy. 
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PACS numbers: 05.40.-a, 91.30.Dk, 87.19. Nn, 05.45.Tp 

On October 17, 2005 two strong earthquakes (EQs) with magnitude around 6.0-units occurred at 05:45:20 and 
12:46:57 UT with epicenters (according to USGS) at 38.15 ! W,26.68 < W and 38.13°7V,26.65°iV, respectively (see FigHJ. 
. ■ At the same epicenter, a third almost equally strong earthquake occurred at 21:40 UT on October 20,2005. All the 
three epicenters lie at a distance of around 100km from Lesvos island, at which the MYT station -on the dipoles of 

• i-H , which the intense signals Mi to M4 (Fig. 1(a), (b) of the main text) have been recorded- is located. This verifes that 
these signals are actually Seismic Electric Signals (SES) as classified in advance (i.e., upon the initial submission of 
the paper on April 16,2005) in the main text on the basis of the entropy criterion proposed. 

We now follow the procedure described in Refs.0,0], in order to investigate whether the time window of the impend- 

L - - '! ing strong EQs could have been estimated. We consider either the region A-.N^qE^^ or the region B-.N^^E^'q, 
which surround the EQ epicenters and the MYT station (see FigUJ, and study how the seismicity evolved after the 
SES initiation. If we set the natural time for seismicity zero at the initiation of the concerned SES activities, we form 

OO . time series of seismic events in natural time for various time windows as the number N of consecutive (small) EQs 

£■»-. ' increases. We now compute the normalized power spectrum 0,0 in natural time H(<f>) for each of the time windows 
and the results are depicted in FigEl As examples we consider in this figure two magnitude thresholds (herafter 
referring to the local magnitude Ml or the 'duration' magnitude Mr>) 3.4 (upper) and 3.6 (lower). In the same figure, 
we plot in blue the power spectrum obeying the relation 

_, . 18 6 cos uj 12 sin u 

U{UJ) =5&- 57^ W 

O . 

which holdsP,0,3 when the system enters into the critical stage {u> — 2n(j), where <f> stands for the natural frequency^, 
£>V 0,0)- An inspection of Fig[21reveals that the red line approaches the blue line as N increases and a coincidence occurs 
at the last small event which had a magnitude 3.6 and occurred at 04:31 UT on October 17, 2005, i.e., roughly one 
hour before the first strong EQ. To ensure that this coincidence is a true oneP,0,@we also calculate the evolution of 
K\,S and S- (cf. Ki stands for the variance K\ = (\ 2 ) — (x) 2 as explained in RefsTP,0]) and the results are depicted 
in Fig|3|for three magnitude thresholds 3.4, 3.5 and 3.6. 

We now further comment on the aforementioned results. Since the strong EQs occurred in the border between Greece 
and Turkey, the seismicity catalogues of neither Greek nor the Turkish Institutes can be considered as complete for 
small magnitudes. Hence, we preferred here to make the calculations by using the United States Geological Survey 
(USGS) catalogue (see Table[IJ|. Irrespective if we use the seismicity in the region N^qE^ ® or in the smaller region 
-/Vly l-Eff q , the coincidence occurs upon the occurence of the aforementioned 3.6 EQ (almost 1 hour before the first 
strong EQ). The magnitude of this EQ comes from the European-Mediterranean Seismological Centre (EMSC) (see 
the corresponding announcement in Fig^J since it has not yet been reported by USGS. Note that if we take the 
magnitude of this EQ to be somewhat larger, then the first box in Fig|3] (which has been plotted for magnitude 
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FIG. 1: Map of the area surrounding the measuring station of MYT and the epicenters of the three strong EQs that occurred 
on October 17 and October 20,2005. The earthquake mechanisms of all three EQs are also shown. The seismicity subsequent 
to the SES initiation has been studied in the gray shaded regions. 
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FIG. 2: The normalized power spectrum(red) H(4>) of the seismicity in area A as it evolves event by event after the initiation 
of the SES activities Mi to M4. The two examples presented correspond to the two different magnitude thresholds 3.4 and 3.6 
in (a),(b) respectively. In each case only the region 4> £ [0,0.5] is depicted (for the reasons discussed in Refs.jlJ, |2J), whereas 
the 11(0) of Eq.Q is depicted by blue color. 



threshold 3.4) shows that the coincidence occurs on the last but one event, i.e., on October 13,2005 (almost three 
days before the first strong EQ). 

The following comment might be useful: In the analysis of signals depicted in Figsl(a),(b) we proceeded as follows: 
We first analyzed M\ (recorded on March 21, see Figl(a)) and found that it obeys the criterion (i.e., S and S*_ smaller 
than S u )■ We then turned to the recordings on March 23: we first considered M4 -which is well distinguishable from 
the others in view of its opposite polarity- and found again that the criterion is fullfillcd. As for the remaining 
recordings of March 23 (comprising M2 and M3) we checked both possibilities that is: (1) we considered M2 and 
M3 together and found that not only the criterion was violated but also the M2 + -M3 signal (altogether) behaved 
like signals obtained from a "uniform" distribution(i.e., K\ was around 1/12). (2)on the other hand, if we analyse 
separately M2 and M3 then both results obeyed the criterion, (thus being consistent with the conclusions drawn from 
Mi and M4). Hence, we preferred to draw in Fig. lb the latter possibility. It has been earlier suggested (see pages 
18 and 114 of Ref.[10] of the main text) that the change of the SES polarity might be associated with a different 
EQ source mechanism. The first three SES activities Mi, M2 and M3 -which have been plotted in Figs. l(a,b) do 
have the same polarity and hence- might be associated with the 3 EQs that have already occurred having the same 
mechanism. If this suggestion is correct, the last signal, i.e., M4, with opposite polarity, should correspond to another 
EQ with different mechanism which has not occurred yet. 



[1] P. A. Varotsos, N. V. Sarlis, and E. S. Skordas, Practica of Athens Academy 76, 294 (2001). 
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FIG. 3: Evolution of Ki, S and S- upon using the USGS catalogue for various magnitude (Ml or Md) thresholds for the two 
regions A and B. 
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TABLE I: The United States Geological Survey (USGS) Catalogue for the area A under discussion together with the last event 
reported by EMSC (see FigEJ. 



Number 


Year 


Month 


Day 


Hour 


Min 


Sec 


Lat 


Lon 


d(km) 


M 


Source 


Ma 


Source 


1 


2005 


3 


30 


13 


57 


46.38 


38.782 


25. 


.782 


56 


3.4 


MDISK 


3.4 


MDATH 


2 


2005 


4 


1 


21 


41 


39.5 


38.73 


25. 


,84 


26 


3.5 


MDATH 






3 


2005 


5 


7 


19 


4 


40.31 


39.492 


26. 


.158 


12 


3.5 


MDATH 


3.5 


MDISK 


4 


2005 


5 


11 


22 


3 


3.5 


38.77 


25. 


,56 


31 


3.4 


MLATH 






5 


2005 


5 


13 


22 


30 


29.7 


38 


26. 


,68 


12 


3.7 


MLTHE 


3.7 


MDATH 


6 


2005 


6 


3 


13 


48 


56.8 


38.78 


26. 


,54 


54 


3.5 


MDATH 






7 


2005 


6 


12 


12 


43 


46.1 


38.01 


27 


,61 


5 


3.5 


MDISK 


3.5 


MDATH 


8 


2005 


6 


17 


2 


47 


54.6 


38.06 


27 


,4 


27 


3.8 


MDATH 






9 


2005 


6 


23 


11 


55 


39.2 


37.8 


26. 


,76 


25 


3.9 


MLATH 






10 


2005 


6 


23 


13 


7 


6.2 


37.78 


26. 


,81 


17 


3.6 


MDATH 






11 


2005 


6 


23 


22 


44 


16.7 


37.76 


26. 


,7 


19 


4.1 


MLTHE 


4.1 


MLATH 


12 


2005 


6 


24 





21 


33 


37.36 


26. 


,87 


22 


3.4 


MDATH 






13 


2005 


6 


24 


4 


27 


8.6 


37.76 


26. 


,72 


45 


3.8 


MDATH 






14 


2005 


6 


26 


20 


7 


26.7 


37.77 


26. 


,77 


18 


3.8 


MLATH 






15 


2005 


6 


27 


17 


59 


51.6 


38.87 


27 


,67 


23 


3.6 


MDATH 


3.4 


MDISK 


16 


2005 


6 


30 


1 


3 


31 


37.68 


26. 


,8 


22 


3.5 


MDATH 






17 


2005 


6 


30 


21 


34 


53.8 


37.8 


26. 


,78 


19 


3.8 


MLATH 






18 


2005 


6 


30 


22 


14 


52.6 


37.73 


26. 


,86 


20 


3.4 


MDATH 






19 


2005 


7 


11 


3 


23 


30.2 


37.79 


27 


,47 


23 


3.7 


MDATH 


3.6 


MDISK 


20 


2005 


7 


12 





47 


45.4 


37.8 


27 


,59 


24 


3.6 


MDATH 


3.3 


MDISK 


21 


2005 


7 


12 





55 


19.5 


37.83 


27 


,6 


10 


3.5 


MDATH 


3.2 


MDISK 


22 


2005 


7 


17 


11 


6 


39.4 


37.75 


26. 


,69 


5 


3.5 


MDATH 






23 


2005 


7 


19 


3 


56 


52 


38.754 


27 


,829 


21 


3.5 


MDATH 


3.4 


MDISK 


24 


2005 


7 


21 


3 


12 


35 


37.774 


27 


,574 


6 


3.4 


MDISK 


3.4 


MDATH 


25 


2005 


7 


25 


21 


56 


48 


38.764 


27 


,721 


17 


3.4 


MDATH 


3.2 


MDISK 


26 


2005 


8 


3 


23 


52 


59 


39.32 


26 


,88 


28 


4.1 


MLATH 


4.1 


MLISK 


27 


2005 


8 


6 


16 


40 


7.7 


38.88 


26. 


,34 


24 


3.4 


MLISK 


3.4 


MDATH 


28 


2005 


8 


9 


9 


11 


11.4 


37.06 


27 


,89 


54 


3.4 


MDATH 


3.2 


MDISK 


29 


2005 


8 


15 


20 


51 


22 


38.084 


26. 


,314 


12 


3.6 


MLATH 


3.5 


MLISK 


30 


2005 


8 


20 


11 


53 


17 


38.599 


26. 


,041 


29 


3.5 


MDATH 


3.3 


MLISK 


31 


2005 


8 


25 


13 


9 


10 


37.222 


28 




24 


3.5 


MDATH 


3 


MDISK 


32 


2005 


8 


30 


17 


39 


46.43 


39.141 


27 


,9 


16 


4 


MLISK 


4 


MDATH 


33 


2005 


9 


1 


19 


24 


56.4 


37.77 


26. 


,75 


19 


3.6 


MDATH 






34 


2005 


9 


4 


18 


14 


37.8 


38.75 


25. 


,59 


38 


3.5 


MDATH 






35 


2005 


9 


11 


23 


30 


47.3 


39.24 


26. 


,62 


24 


3.4 


MDATH 






36 


2005 


10 


13 


15 


33 


8.2 


38.32 


26. 


,84 


25 


3.4 


NOA 






37 


2005 


10 


17 


4 


31 


27.3 


38.19 


26. 


,71 


20 


3.6 


EMSC 
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EMSC relocation in GSE format 
BEGIX CSB2.0 
MSG_{}TYPE DATA 
MSG_{}ID 051017043127 EMSC 
DATA {}TYPE BULLETIN GSE2.0 
EVENT 20051017043127 

Date Time Lath ude Longitude Depth Ndef Nsta Gap Magi N Mag2 N Mag3 N Author ID 
rms OT_{}Error Smajor Sminor Az Err mdist Mdist Err Err Err Quality 
2005/10/17 04:31:27.3 38.1899 26.7128 20.0 44 41 80 ML 3.6 CSEM GSP66 
0.76 +- 0.12 3.4 2.5 84 0.32 5.90 a i uk 
AEGEAN SEA 

Sta Dist EvAz Phase Date Time TRes Azim AzRes Slow SRes Def SNR Amp Per Magi Mag2 ID 

BLCB 0.32 52.9 m Pg 2005/10/17 04:31:35.1 0.3 T KAN 

BLCB 0.32 52.9 m Sg 2005/10/17 04:31:40.1 -0.2 T KAN 

SMG 0.49 168.4 a Pn 2005/10/17 04:31:36.6 -0.1 T NOA 

IZM 0.49 64.7 m Pg 2005/10/17 04:31:37.8 -0.1 T KAN 

AKS 1.10 51.1 m Pn 2005/10/17 04:31:48.2 -0.1 T KAN 

PRK 1.11 342.0 a Pn 2005/10/17 04:31:47.9 -0.1 T NOA 

BODT 1.22 157.0 m Pn 2005/10/17 04:31:49.6 -0.3 T KAN 

BODT 1.22 157.0 m Sg 2005/10/17 04:32:06.6 -2.2 T KAN 

MLSB 1.23 136.3 m Pn 2005/10/17 04:31:50.7 0.6 T KAN 

APE 1.46 220.4 m Pn 2005/10/17 04:31:51.3 -1.6 T KAN 

YER 1.63 129.6 m Pn 2005/10/17 04:31:55.7 -0.2 T KAN 

EZN 1.66 349.7 m Pn 2005/10/17 04:31:56.0 -0.2 T KAN 

BALB 1.71 31.8 m Pn 2005/10/17 04:31:56.7 -0.2 T KAN 

CANB 1.84 8.4 m Pn 2005/10/17 04:31:57.7 -1.1 T KAN 

DENT 1.88 102.6 m Pn 2005/10/17 04:31:59.8 0.4 T KAN 

DST 2.05 45.9 m Pn 2005/10/17 04:32:01.7 -0.1 T KAN 

LIA 2.08 325.5 a Pn 2005/10/17 04:32:02.0 0.6 T NOA 

DALT 2.09 132.2 m Pn 2005/10/17 04:32:02.5 0.2 T KAN 

LPK 2.18 1.0 m Pn 2005/10/17 04:32:03.8 0.2 T KAN 

KHL 2.22 85.7 in Pn 2005/10/17 04:32:04.9 0.6 T KAN 

EDO 2.33 22.2 m Pn 2005/10/17 04:32:06.1 0.4 T KAN 

BNT 2.35 22.6 m Pn 2005/10/17 04:32:05.0 -1.0 T KAN 

FETY 2.44 128.7 m Pn 2005/10/17 04:32:07.2 2.8 T KAN 

ORLT 2.51 41.8 m Pn 2005/10/17 04:32:08.4 0.0 T KAN 

MRMT 2.51 15.6 m Pn 2005/10/17 04:32:08.6 0.3 T KAN 

MFT 2.63 9.4 m Pn 2005/10/17 04:32:09.3 -0.8 T KAN 

KARP 2.66 1 72.1 a Pn 2005/10/17 04:32:09.3 -0.2 T NOA 

ALT 2.80 71.0 in Pn 2005/10/17 04:32:12.3 -0.2 T KAN 

ELL 2.92 118.5 m Pn 2005/10/17 04:32:14.7 0.4 T KAN 

ISP 3.03 95.8 mdi P 2005/10/17 04:32:15.3 -0.3 T GSRC 

ISP 3.03 95.8 m Pn 2005/10/17 04:32:16.0 0-2 T KAN 

NPS 3.05 197.2 a Pn 2005/10/17 04:32:15.6 0.8 T NOA 

RDO 3.09 343.3 a Pn 2005/10/17 04:32:15.8 0.5 T NOA 

YLV 3.14 40.1 in Pn 2005/10/17 04:32:17.8 0.4 T KAN 

BOK 3.15 102.2 m Pn 2005/10/17 04:32:16.8 -0.8 T KAN 

CTT 3.24 23.6 m Pn 2005/10/17 04:32:18.4 -0.3 T KAN 

BADT 3.24 34.2 m Pn 2005/10/17 04:32:18.4 -0.3 T KAN 

VLI 3.34 245.1 a Pn 2005/10/17 04:32:20.0 1-2 T NOA 

ESKT 3.49 66.3 m Pn 2005/10/17 04:32:22.7 0.2 T KAN 

EDRB 3.65 0.4 m Pn 2005/10/17 04:32:24.5 -0.1 T KAN 

KIZT 4.11 78.7 m Pn 2005/10/17 04:32:31.2 -0.1 T KAN 

MDI] 4.17 55.5 m Pn 2005/10/17 04:32:32.4 0.6 T KAN 

BRTR 5.61 72.0 m i P 2005/10/17 04:32:52.8 0.6 T GSRC 

IKL 5.90 107.2 m Pn 2005/10/17 04:32:56.6 0.1 T KAN 

STOP 



FIG. 4: The detailed European-Mediterranean Seismological Centre (EMSC) report for the 3.6 EQ that occured almost one 
hour before the first strong EQ. 



